Harmonic entanglement in a degenerate parametric down conversion 
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We study the harmonic entanglement and squeezing in a two-mode radiation produced in a degen- 
erate parametric down conversion process coupled to a two-mode vacuum reservoir employing the 
linearization procedure. It is found that there is a quadrature entanglement between the harmon- 
ically related fundamental and second-harmonic modes and the superimposed radiation exhibits a 
significant two-mode squeezing. The entanglement exits even when there is no two-mode squeezing, 
since the correlation leading to these phenomena are essentially different. In addition, the more the 
external coherent light is down converted, the more stronger the entanglement and mean photon 
number of the two-mode radiation would be which is not generally true for squeezing. 

PACS numbers: 42.50.Dv, 42.65.Yj, 03.65.Ud 



I. INTRODUCTION 

Optical degenerate parametric down convrsion is one 
of the second-order nonlinear processes in which a pump 
photon of frequency 2lu is down converted into a pair of 
signal photons each of frequency u>. Due to the inher- 
ent two-photon nature of the interaction, the paramet- 
ric oscillator or amplifier is one of the most interesting 
and well studied devices in the nonlinear quantum optics 
EE S S, i, 1, !, 0, ©, 0- It is found to be a good 
source of light with nonclassical features. The quantum 
optical properties are significantly degraded by the leak- 
age through the mirrors and amplification of the quan- 
tum fluctuations in the cavity. However, if the ordinary 
vacuum is replaced by a squeezed vacuum reservoir, the 
nonclassical properties can be enhanced pro- 
vided that the reservoir is squeezed in the right quadra- 
ture, since the fluctuations entering the cavity is biased. 
Moreover, when a nonlinear crystal is shined with an ex- 
ternal radiation of frequency 2u> only some part of this 
radiation is down converted into a pair of photons Q. 
As a result, the cavity contains the down converted and 
unchanged radiations. Although a considerable attention 
has been given to a single-mode squeezed radiation pre- 
viously, most recently it is envisioned that a degenerate 
parametric oscillator can be a source of a two-mode light 
with nonclassical properties characterized by a strong 
correlation between the states of each mode separately 
as well as the superimposed state formed by the two ra- 
diations 

Currently, generation of macroscopic entangled states 
is receiving attention in connection with their poten- 
tial in the quantum information and measurement the- 
ories [ill EH- At microscopic level nonclassical corre- 
lations and entanglement arise due to the down con- 
version of a single high frequency photon into a pair 
of correlated lower frequency photons. The generation 
of the entangled beams in the nondegenerate paramet- 



*Electronic address: sint'tesfa@yahoo.com 



ric down conversion was predicted by Reid and Drum- 
mond 0, HH, HE]- Studies have also shown the exis- 
tence of the correlation between the fundamental and 
second-harmonic modes [TtJ including entanglement in 
the second- harmonic generation (l8j . In this respect, Lim 
and Saffman [19j analyzed the production of two beams 
with nonclassical intensity correlations and quadrature 
entanglement in the dual-ported reservoir of the second- 
harmonic generation and found that the harmonic output 
exhibits strong quantum correlations. In connection to 
this, Grosse et al. [2(| predicted a perfect entanglement 
between the fundamental and second-harmonic modes in 
the pump depleted nondegenerate parametric amplifica- 
tion and defined such an entanglement between harmon- 
ically related fields as harmonic entanglement. 

Einstein-Podolsky-Rosen (EPR) presented in 1935 [U 
their famous argument, which on the basis of the local 
realism claimed that an observation of perfectly corre- 
lated positions and momenta would imply the incom- 
pleteness of the quantum theory. The essential step in 
realizing the EPR-type entanglement is then to intro- 
duce correlated states of at least two particles that per- 
sist even when the particles are spatially separated. As 
thoroughly discussed by Dechoum et al. [22[ without 
having to construct an experimentally impossible states 
with perfect correlation, as in the original paper, it is pos- 
sible to demonstrate this type of correlation using the 
inferred Heisenberg principle. In this regard, a direct 
and experimentally feasible qualitative criterion for such 
a correlation of continuous variables was first proposed 
by Reid [l5| via quadrature phase amplitude, applying a 
nondegenerate parametric amplification, which is closely 
related to the original version. In addition, making use 
of the quadrature variables Lodahl [23j showed the exis- 
tence of EPR-type correlations in the second-harmonic 
generation. Most recently, Olsen has considered 
a travelling-wave second-harmonic generation and pre- 
dicted that this quantum system can be employed in 
the experimental demonstration of entanglement with 
continuous variables. It then appears natural to ask 
whether the strong correlation between the fundamental 
and second-harmonic modes of the degenerate parametric 
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oscillator associated with the down conversion also leads 
to entanglement or not? The main task of this commu- 
nication is, therefore, devoted to answer this question. 

In this paper, the nonclassical properties of the two- 
mode cavity radiation of the driven degenerate paramet- 
ric oscillator coupled to a two-mode vacuum reservoir 
applying the linearization procedure in which the quan- 
tum properties of the system in time are taken to vary 
slightly around the steady state mean values would be 
analyzed. As discussed by Gilles et al. this ap- 

proximation remains valid as long as the quantum fluc- 
tuations are much smaller than the classical mean values 
which is the essence of a weaker coupling of the radiation 
with the oscillator that corresponds to small quadrature 
noise. We also consider the semiclassical approximation 
whereby the two modes are assumed to be uncorrelated 
at steady state. As recently reported by Chaturvedi et 
al. [25l | , the semiclassical theory is found to work surpris- 
ingly well in the threshold region. We hence study the 
harmonic entanglement and its relation with the squeez- 
ing near threshold, since the relation between the two 
is an interesting issue by its own right [26|, [27[. More- 
over, we investigate the efficiency with which a nonlinear 
crystal down convert the light falling on it and the associ- 
ation of the down conversion with the entanglement and 
squeezing. We also calculate the mean photon number 
to see how intense the generated light could be. 



II. QUANTUM LANGEVIN EQUATIONS 

Interaction of an external coherent radiation with a 
nonlinear crystal responsible for a degenerate parametric 
oscillation placed in a resonant cavity can be described 
in the rotating-wave approximation and in interaction 
picture by the Hamiltonian of the form 



is 



P - b 



(1) 



where e is proportional to the amplitude of the coher- 
ent input, A is the measure of the coupling of a nonlin- 
ear crystal with the external coherent radiation, a and 
b are the time-independent annihilation operators for 
the fundamental (subharmonic) and pumping (second- 
harmonic) modes, and A and e are chosen to be real- 
positive constants. In view of the fact that a and b are 
mutually commuting operators, the pertinent quantum 
Langevin equations are found to be 



- = Xa^--a + F a (t), 



(2) 



Langevin noise operators satisfying, for a two-mode vac- 
uum reservoir, the correlation functions: 



(£(*)>=«, 



(4) 



(^(t)^(f )> = (FRt)FUt')) = (F^Ftf)) = 0, (5) 



(Fi(t)Fkt'))i^=0- 



(6) 



(7) 



We notice that Eqs. ([2]) and ([3]) are nonlinear coupled 
differential equations that can be solved employing the 
linearization procedure. In this approach, we first take 



a(t) = a + A(t), 



b(t) = p + B(t), 



(8) 



(9) 



where A(t) and B(t) are very small variations about the 
mean values at steady state, and we define a = (a(t)) ss 
and B — (b(t)) ss . This approximation remains valid as 
long as the quantum fluctuations about the mean values 
are much smaller than the classical mean values, which 
implies that the quantum noise during the interaction is 
quite small. Upon taking the statistical average of Eqs. 
([2]) and ([3]) and then using the semiclassical approxima- 
tion, whereby at steady state the modes are assumed 
to be uncorrelated, (a' (t)b(t)) ss — (a'(t)) ss (b(t)) ss , and 
classical decorrelation that operators are assumed to be 
factorized, (a 2 (t)) ss = (a(t))l s , we obtain 



Xa* a 

2 



Aa 2 + Kj3 = 2s. 



(10) 



(11) 



We realize that the semiclassical assumption is found to 
work for weak nonlinearity or weak coupling between the 
external radiation and nonlinear crystal where the mean 
photon number at threshold is very large. 

Now with the aid of Eqs. ©, ©, ©, ®, (HJJJ, (D), 
and the fact that A and B are small perturbations, it is 
possible to verify that 



dA(t) 
dt 



: e\B(t) + e 2 A\t) 



%A(t)+F a {t), (12) 



dt 



—a 

2 



F b (t), 



(3) 



where k is the cavity damping constant chosen to be the 
same for both modes and Fi(t), with i = a, b, are the 



dB(t) 
dt 



- Sl A(t)- -B(t)+F b (t), 



(13) 



where we set = Aa* and £2 — A/3. Multiplication of 



Eq. (jlO)) by a suggests that the phase of 3 is the same 
as that of a 2 . This is consistent with Eq. (fTTj) if both 
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a and f3 are real, which implies that e\ and £2 are also 
real. It then follows from Eqs. (flU|) and (fTTj) that 



(14) 



and £2 = ft/2 for £1 7^ 0. We notice that there are two 
possible values for the fundamental amplitude and hence 
the system may be in a transient su per imposed state of 
these amplitudes prior to detection [13|]. Moreover, the 
solutions of these coupled differential equations are found 
following a straight forward algebra to be 

A(t) = oi(t) + [a 3 (t) +pa 4 (t)]A(Q) + [a 5 (t) + pa 6 (t)] 
x it(0) + qa 6 {t)B(0) + qa A {t)&{Q) + f(t), (15) 

B{t) = a 2 (t) + [03 (t) - pa A {t)} B(0) + [a 5 (t) - pa 6 (t)] 
x flt(0) - qa 6 (t)A(Q) - ga 4 (t)i*(0) + (16) 

where 



/(f) = / [(a 3 (t-0+Pa4(t-0)^a(*') 
Jo 

+ (a 5 (t-t')+pa 6 (t-t'))F^t') 

+ qa 6 (t - t')F b (t') + qa A (t - t')F^{t')]dt' , (17) 



[(a 3 (t~t')-pa 4 (t-t'))F b (t') 

+ (a 5 {t-t')- P a 6 (t-t'))F b i (t') 

- qa 6 (t - t')F a (t') - qa 4 (t - t')F^t')]dt' , (18) 



in which 



£2 



2ei 



(19) 
(20) 



Ol (t) = — 



£1 cosh 



*V^4£l 



(£ip + £ 2 g) sinh 



(21) 



(*-E 2 )« 



02(t) 



(eig + £2^) sinh 



- 4 E ? 



-£2 cosh 



ty/e'i - 4e( 



(22) 



«:<(/) = r - cosh (^) COsh M j , ( 23 ) 



a^it) — e ~ 2 " sinh ( ^ 



a5(t) = e ^ sinh ( ^ 



a 6 (t) = e cosh ( 



sinh 



cosh 



sinh 



4£? 



- 4£^ 



(24) 



(25) 



(26) 



It may worth to mention that Eqs. (|T5|) and (|T6|) are 
applied to calculate various quantities of interest. We 
also observe that the condition for e\ to be real requires 
that 2Ae > K£2- We, therefore, denote the case for which 
2Ae = «£ 2 that corresponds to £1 = as a threshold 
condition. 



III. QUADRATURE ENTANGLEMENT 

In this section, we seek to study the entanglement of 
the fundamental and second-harmonic modes in the cav- 
ity. It is a well-established fact that a quantum system 
is said to be entangled, if it is not separable. That is, 
if the density operator for the combined state cannot be 
expressed as a combination of the product density oper- 
ators of the constituents, 



(27) 



On the other hand, entangled continuous variable state 
can be expressed as a co-eigenstate of a pair of EPR-type 
operators such as X a — Xb and P a + P b [28j]. The total 
variance of these two operators reduces to zero for maxi- 
mally entangled continuous variable states. Nonetheless, 
according to the criterion set by Duan et al. [291 ] quan- 
tum states of the system are entangled, provided that the 
sum of the variances of a pair of EPR-likc operators 



u = X a - X b 



v=P,.+ Ph. 



(28) 



(29) 



where X a = -± (6 + fit) , X b = -± (b + 6+) , 

^ a = 72 - ®) ' and ^ b = 72 - ^) ' satisf y 



Au 2 + Av 2 < 2, 



(30) 



in which 



Au 2 + Av 2 = 2(6*6) + 2(6*6) + 2(a)(6) + 2(6*) (6*) 

- 2(66) - 2(6*6*) -2(6*) (6) - 2(6*) (6) + 2. 

(31) 
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Next we determine the various correlations in Eq. (|3Tj) . 
To this end, assuming the cavity mode to be initially in 
the vacuum state and using the fact that the Langevin 
noise forces have zero mean along with Eqs. ||SJ), ©, 
CPU), and (HU), we get 



(a(t)) =a + oi(t), 



(S(t))=/3 + a 2 (t), 



(32) 



(33) 



(at(i)s(t)) = + «iW) 2 + (/ f (*)/(*)>, (34) 



(6t(i)S(i)) = (/3 + a2 (t)) 2 + (gt (W))i (35) 



(o(*)6(t)) = (a + ai (t))(/3 + a 2 (t)) + (f(t)g(t)). (36) 

Moreover, in view of the correlations of the Langevin 
noise forces (O, ([5]), and ([7]), it is possible to show at 
steady state that 



</*(*)/(*)>« = 



k[(1 + p 2 + g 2 )(/t - eg) - 2 Pv / £ 2 _ 4e 2 ] 
8[«(«- 2e 2 ) +4e 2 ] 



k[(1 +p 2 + <? 2 )(k + eg) + 2 Pv /^ - 4e 2 ] 
8[k(k + 2e 2 ) + 4ef] 

K 2 (l -p 2 - (J 2 ) K 2 (l-p 2 +gr 2 ) 



4(k 2 - £ 2 ) 4[k 2 - ( £ 2 - 4s 2 )] 
(l+p 2 -q 2 ) 



(37) 



k[(1 +p 2 + g 2 ){ K - eg) + 2p^£ 2 - 4e 2 ] 
8[k(k-2 £2 )+4£ 2 ] 



k[(1 +p 2 + 9 2 )(k + e 2 ) - 2p^e 2 - 4e 2 ] 

8[k(k + 2 £2 ) +4e 2 ] 
K 2 (l- p 2 - q 2 ) K 2 (l-p 2 + q 2 ) 

4(k 2 - 4) 
(l+p 2 -q 2 ) 



4 [ K 2 _ (e 2 - As 2 ) 



(38) 



(f(t)g(t)) s 



npq(n - e 2 ) 
~4[k(k-2 £2 ) + 4£ 2 ] 

k<jp(k + e 2 ) 
4[k(k + 2£ 2 ) + 4£ 2 ] + 2(k 2 -e 2 ) 



qpn 



K 2 pq 



pq 



(39) 



2[k 2 -( £ 2 -4 £ 2 )] 2" 
We, therefore, see with the aid of Eqs. (|5T]). {35]), ([3"3")) . 



ill, dsa 

Au 2 



36|). (|37 |) . (|38 j) . and (|39 [) at steady state that 
. , n(l+p 2 + q 2 + 2qp){ K ~e 2 ) 



Av 



2[k(k- 2£ 2 ) + 4e 2 ] 
K(l+p 2 +g 2 + 2gp)(K + £ 2 ) 
2[k(k + 2£ 2 ) + 4£ 2 ] 



K 2 (l -p 2 



2qp) 



(k 2 - £ 2 ) 

K 2 (l -p 2 + q 2 - 2qp) 

[« 2 -(£ 2 -4£ 2 )] 

(l+p 2 -g 2 + 2^). 



(40) 



Finally, on account of Eqs. (flT))) . (|2H)l . and the fact that 
£ 2 = re/2 for £i 7^ 0, we find 



Aw 2 + Av 2 



+ 



k 3 (k + 4ei) 
8e 2 (k 2 - 16e 2 ) 

3k 3 (k + 4ei) 
4(k 2 -16£ 2 )(k 2 + 2£ 2 ) 
16ei(4ei + re) + 24k£i 
3(k 2 - 16£ 2 ) 
16k 3 £i 
(k 2 -16£ 2 )(3k 2 + 16£ 2 ) 



(41) 



It is not difficult to observe that Au 2 + Av 2 is very large 
when £i = and E\ = 0.25k. We hence plot Am 2 + Av 2 
versus £i for ei > 0.25k. 



K=0.1 

k=0.2 
k=0.3 




FIG. 1: Plots of the sum of the variances of the EPR-type 
operators (Au 2 + Av 2 ) at steady state for different values of 



As clearly shown in Fig. 1 the correlation between 
the fundamental and second-harmonic modes exhibits a 
quadrature entanglement except near certain value of e\, 
for example e\ — 0.13 for k = 0.3 at steady state. We also 
realize that entanglement does not exist near a thresh- 
old value, £i = 0. However, entanglement exists and 
it decreases with damping constant in other cases for 
which £i > 0.25k. This is related to a well known fact 
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that the lesser the cavity damping constant, the more 
the radiation stays in the cavity which in turn enhances 
the correlation that leads to entanglement. It is possible 
to observe that the dependence of the entanglement on 
damping through the mirrors is insignificant for larger 
values of £\. Moreover, as can easily be seen from Eq. 
(pTTj) the existing entanglement is attributed to the cor- 
relation between the states of the two cavity modes. We 
realize that since E\ — X(a) ss , it corresponds to the de- 
gree at which the external radiation of frequency 2to is 
down converted by the nonlinear crystal. Therefore, it 
can be inferred from Fig. 1 that the entanglement would 
be stronger, the more efficiently the external radiation is 
down converted by the crystal. This indicates that even 
though the down conversion process breaks the coherent 
external radiation into two, it is unable to destroy the 
coherence that is responsible for the correlation between 
the down converted and unchanged radiations. 
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FIG. 2: Plots of the minus quadrature variance for the two- 
mode cavity radiation (Ac?.) at steady state for different val- 
ues of K. 



IV. QUADRATURE VARIANCES 



In this section, we seek to analyze the squeezing prop- 
erties of the two-mode cavity radiation which can be de- 
scribed by annihilation operator 



(42) 



where fi and b are the boson operators that represent 
the fundamental and second- harmonic modes. In view 
of the boson commutation relation for d and b, one can 
easily see that [c, cH = 1 and [c, c] = 0. We note that 
the squeezing of the two-mode cavity radiation can be 
studied using the quadrature operators corresponding to 



c+ = cr + c 



and 



(ct-c), 



(43) 
(44) 



in which the squeezing occurs when 

Ac 2 ± = 1 + (d f d) + (tfb) ± ((ab) + (a+St)) + (a+fo) 

+ <aS+) ± |[(a+ 2 ) + (5t 2 > + (a 2 > + (6 2 ) + (a) 2 

+ (S) 2 + <at) 2 + (St) 2 ] ± ((d)(6) + <at><St» 
+ <o><a+> + <a><St> + <S><at> + <S><8t> (45) 



is less than one. Following a similar approach as in Sec- 
tion III, it is possible to see that 



k 2 (k =f 2e 2 ) + nel + Ake\ ± 2k£ 2 £i 



Ac z ± = 



(kT£ 2 )[k(kT 2e 2 ) +4e 2 ] 



(46) 



which can also be put on the basis of the fact that e 2 
k/2 for Ei ^ in the form 



K 2( 5=F 4)_|_4 ei (4 £l±K ) 



AC ± 2(2 T l)(« 2 (lTl)+4£ 2 ) 



(47) 



It turns out that the superimposed cavity radiation ex- 
hibits a significant two-mode squeezing for certain values 
of Ei corresponding to each damping constant. As clearly 
presented in Fig. 2 the degree of squeezing decreases with 
the cavity damping constant for smaller values of Si, but 
it increases for larger values. It is not difficult to see that 
a maximum obtainable squeezing occurs slightly above 
the threshold value and it is independent of the cavity 
damping constant. It is also found that a two-mode max- 
imum squeezing of about 42% occurs at different values 
of Ei for different damping constants. The value of e± for 
which a maximum squeezing occurs increases with the 
cavity damping constant. Comparing the results given 
in Figs. 1 and 2 shows that two-mode squeezing exits 
even for values of e± for which there is no entanglement. 
It goes without saying that the differences in the entan- 
glement and squeezing properties are essentially related 
to the correlations that lead to these phenomena. That 
is, the correlations between similar states of a radiation 
like (d 2 ) and (6 2 ) contribute to the squeezing, but not to 
the entanglement, since entanglement requires two differ- 
ent states of radiation to be correlated. 



MEAN OF THE PHOTON NUMBER AND 
INTENSITY DIFFERENCE 



The mean photon number of the cavity radiation cor- 
responding to the superposition of the two modes can be 
expressed as 



(t)c(t)), 



(48) 
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where c(t) is the annihilation operator defined by Eq. 
([4"2"]l . Hence it is not difficult to see that 



R = -[<at(t)a(t)> + <St(*)6(t)) 
+ <at(*)S(t)> + <S+(t)a(*))]. 



(49) 



In order to obtain the mean photon number, we need to 
determine the involved cross correlations. To this effect, 
making use of Eqs. ®, ®, ([IS]), (TSJ), Q3>, OH), PI- 
(|2T)]) . and the correlations of the Langevin noise forces, 
we find at steady state 



(a* (*)&(*)}„ = (tf(t)a(t))ss = 

K - e 2 



4(ei-4£?) 
k + e 2 2e 2 



k(k-2£ 2 ) + 4e 2 «(« + 2e 2 ) + 4ef 



k 2 - £9 



(50) 



Thus, applying Eqs. 
we arrive at 



, d3SD, (EH), (Hi, and 0, 



1/ | «(«-£ 2 )(l+p 2 +g 2 ) 

2 ^* + /j t 8 K/c-2£ 2 )+4£?] 

k(k + £ 2 )(1 +p 2 + <? 2 ) n 2 (l-p 2 -q 2 ) 



8[k(k + 2s 2 ) + 4e 2 } 



4(k 2 



K 2 (l-p 2 + g 2 ) (l+p 2 -? 2 ) 



K£i£ 2 



4[k 2 



-(£ 2 -4£ 2 ) 
K £ 2 



4 

K + £ 2 



4(£ 2 -4£ 2 ) 

2£ 2 



- 2e 2 ) + 4£ 2 k(k + 2e 2 )+4£ 2 



which can also be put for e\ ^ in the form 



(51) 



(2ei 



4A 2 



k 3 (k - 2ei) 



8A 2 32£ 2 (k 2 - 16£ 2 ) 

3k 3 (k + 2ei) 2si(k + 16ei) 



16(k 2 -16£ 2 )(k 2 + 2 £ 2 ) 3(k 2 -16£ 2 ) 



(52) 



We see from Eq. ([52]) that the mean photon number 
takes large values where there is no entanglement. For- 
tunately, the mean photon number increases with s\ for 
larger values of £\. Therefore, as we can see from Fig. 3 
a considerably intense entangled two-mode light can be 
generated from the driven degenerate parametric oscilla- 
tor, which we believed is an encouraging result. 

On the other hand, the intensity difference can be de- 
fined as 



AI 



a) a - 



b r b. 



(53) 



Upon employing Eqs. ([53]). J35J), (37]), and (|55|). the mean 
of the intensity difference at steady state turns out to be 



AI 



2[k(k- 2e 2 ) +4e 2 ] 



4£ 2 



2[k(k + 2 £2 ) +4£ 2 ] 



(54) 




0.5 1 



1.5 2 2.5 3 3.5 



FIG. 3: Plots of the mean of the photon number (n) at steady 
state for A = 0.5 and different values of k. 



in which using Eq. p^]) along with the fact that £ 2 = k/2 
for £ ^ leads to 



AI 



4e 2 - k 2 
AX 2 



16£ 2 8(k 2 + 2e 2 )' 



(55) 




FIG. 4: Plots of the mean of the intensity difference (A7) at 
steady state for A = 0.5 and different values of k. 

According to the results shown in Fig. 4, the mean 
of the intensity difference is nonnegative for all values of 
e%. This indicates that, in principle, a degenerate para- 
metric oscillator can down convert a significant amount 
of light falling on it. We notice that the number of down 
converted photons available in the cavity as compared 
to the number of photons of the driving external coher- 
ent radiation depends on the damping constant, since the 
fundamental and second-harmonic modes are assumed to 
leak through the mirrors in the same way. The mean in- 
tensity difference deceases for smaller values of £i, but it 
increases for larger values. On the basis of Eq. (|14p we 
note that £\ increases with the amplitude of the external 
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coherent radiation provided that the cavity damping and 
coupling constants are taken to be independent of the 
amplitude of the driving radiation. In view of this fact 
we observe that the process of parametric down conver- 
sion decreases with the amplitude of the external coher- 
ent radiation close to the critical point, whereas it will 
be enhanced with the external radiation after a particu- 
lar value of £1 which increases with the cavity damping 
constant. 



VI. CONCLUSION 

We study the squeezing properties of the two-mode 
superimposed radiation and entanglement in the funda- 
mental and second-harmonic modes of a driven degener- 
ate parametric oscillator coupled to a two-mode vacuum 
reservoir. It turns out that the cavity radiation exhibits 
a significantly intense two-mode squeezing under certain 
conditions pertaining to the rate at which the external co- 
herent radiation is down converted. A maximum squeez- 
ing of about 42% occurs slightly above a threshold value 
for different amplitudes of the external driving radiation 
for different damping constants. On account of Eq. (T55|) 
we notice that the mean of the intensity difference would 
be positive, provided that more than 33.3% of the ex- 
ternal coherent radiation is down converted. Therefore, 
upon comparing Figs. 2 and 4 we see that a two-mode 



squeezing exists when greater than 33.3% of the exter- 
nal coherent radiation is down converted and maximum 
squeezing occurs when S\ — n. We come to understand 
that although the squeezing of the subharmonic mode has 
attracted a great deal of attention in previous studies, 
the squeezing of the two-mode radiation is also equally 
significant with additional possible applications. 

The strong correlation between the fundamental and 
second-harmonic modes due to their harmonic relation 
leads not only to quadrature squeezing, but also to en- 
tanglement. Unlike the squeezing, the entanglement in- 
creases with the rate at which the external coherent ra- 
diation is down converted. Since the radiation has more 
chance to oscillate back and forth in the cavity for smaller 
damping constant, the probability that it can be entan- 
gled is better. This must be the reason for the decre- 
ment of the entanglement with the cavity damping con- 
stant, which is also true for the squeezing near threshold. 
Comparison of Figs. 1 and 4 shows that contrary to the 
squeezing the entanglement is found to be better in a 
region where the down conversion of the external coher- 
ent radiation is enhanced. These differences in squeez- 
ing and entanglement are attributed to the differences in 
the correlation leading to these phenomena. It is a well 
established fact that entanglement requires a correlation 
between two different states of the radiation, a restriction 
that does not necessarily apply to harmonically related 
radiations. 
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